* Replication Do file for:
* Fox, Haynes, Zellman (2026) "Governmental Discrimination against Religious Minorities in Sub-Saharan Africa, 1990-2023" Africa Spectrum
* DOI: https://doi.org/10.7910/DVN/CVDSL1

*************************************************************************
* Table 1. Correlates of Government Religious Discrimination, 1990-2023 *
*************************************************************************

asdoc reg mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur emajrelnewtchr emajrelnewtmus mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) replace

asdoc reg mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year if emajrelnewtchr==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year if emajrelnewtmus==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) keep(l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur emajrelnewtchr emajrelnewtmus mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras _cons) stat(F) append

*********************
* Robustness Checks *
*********************

**************************************************************************************
* Table A1. Descriptive Statistics of Dependent and Independent Variables, 1991-2023 *
**************************************************************************************

sum mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur emajrelnewtchr emajrelnewtmus mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras year if !missing(mmxx)&!missing(l.wzsocdisx)&!missing(l.lxx)&!missing(yemajrel_pc_rcs)&!missing(ygroup_pc_rcs)&!missing(ygdp_pcap_log)&!missing(ypop_log)&!missing(yv2x_polyarchy)&!missing(yv2regdur)&!missing(emajrelnewtchr)&!missing(emajrelnewtmus)&!missing(mz00mintchr)&!missing(mz00mintmus)&!missing(mz00mintanimist)&!missing(ypolarize_ras)&!missing(year)

*********************************************************************************************
* Table A2. Pearson’s R Correlation Coefficients for non-dichotomous variables, Full Sample *
*********************************************************************************************

pwcorr mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur ypolarize_ras  if !missing(mmxx)&!missing(l.wzsocdisx)&!missing(l.lxx)&!missing(yemajrel_pc_rcs)&!missing(ygroup_pc_rcs)&!missing(ygdp_pcap_log)&!missing(ypop_log)&!missing(yv2x_polyarchy)&!missing(yv2regdur)&!missing(emajrelnewtchr)&!missing(emajrelnewtmus)&!missing(mz00mintchr)&!missing(mz00mintmus)&!missing(mz00mintanimist)&!missing(ypolarize_ras)&!missing(year), sig

***********************************************************************************************************
* Table A3. Pearson’s R Correlation Coefficients for non-dichotomous variables, Christian Majority Sample *
***********************************************************************************************************

pwcorr mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur ypolarize_ras  if emajrelnewtchr==1&!missing(mmxx)&!missing(l.wzsocdisx)&!missing(l.lxx)&!missing(yemajrel_pc_rcs)&!missing(ygroup_pc_rcs)&!missing(ygdp_pcap_log)&!missing(ypop_log)&!missing(yv2x_polyarchy)&!missing(yv2regdur)&!missing(emajrelnewtchr)&!missing(emajrelnewtmus)&!missing(mz00mintchr)&!missing(mz00mintmus)&!missing(mz00mintanimist)&!missing(ypolarize_ras)&!missing(year), sig

********************************************************************************************************
* Table A4. Pearson’s R Correlation Coefficients for non-dichotomous variables, Muslim Majority Sample *
********************************************************************************************************

pwcorr mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur ypolarize_ras  if emajrelnewtmus==1&!missing(mmxx)&!missing(l.wzsocdisx)&!missing(l.lxx)&!missing(yemajrel_pc_rcs)&!missing(ygroup_pc_rcs)&!missing(ygdp_pcap_log)&!missing(ypop_log)&!missing(yv2x_polyarchy)&!missing(yv2regdur)&!missing(emajrelnewtchr)&!missing(emajrelnewtmus)&!missing(mz00mintchr)&!missing(mz00mintmus)&!missing(mz00mintanimist)&!missing(ypolarize_ras)&!missing(year), sig

*****************************************************************
* Table A5. Variable Inflation Factor for Independent Variables *
*****************************************************************

quietly reg mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur emajrelnewtchr emajrelnewtmus mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year, cluster(yccode_group) robust

vif

quietly reg mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year if emajrelnewtchr==1, cluster(yccode_group) robust

vif

quietly reg mmxx l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year if emajrelnewtmus==1, cluster(yccode_group) robust 

vif

************************************************************************
* Table A6. Alternative Specifications of Models, with only SRD or GRD *
************************************************************************

asdoc reg mmxx l.wzsocdisx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur emajrelnewtchr emajrelnewtmus mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) reset

asdoc reg mmxx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur emajrelnewtchr emajrelnewtmus mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.wzsocdisx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year if emajrelnewtchr==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year if emajrelnewtchr==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.wzsocdisx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year if emajrelnewtmus==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras i.year if emajrelnewtmus==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) keep(l.wzsocdisx l.lxx yemajrel_pc_rcs ygroup_pc_rcs ygdp_pcap_log ypop_log yv2x_polyarchy yv2regdur emajrelnewtchr emajrelnewtmus mz00mintchr mz00mintmus mz00mintanimist ypolarize_ras _cons) stat(F) append

*****************************************************************
* Table A7. Minimal Models including Alternative Specifications *
*****************************************************************

asdoc reg mmxx l.wzsocdisx l.lxx i.year, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) reset

asdoc reg mmxx l.wzsocdisx i.year, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.lxx i.year, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.wzsocdisx l.lxx i.year if emajrelnewtchr==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.wzsocdisx i.year if emajrelnewtchr==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.lxx i.year if emajrelnewtchr==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.wzsocdisx l.lxx i.year if emajrelnewtmus==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.wzsocdisx i.year if emajrelnewtmus==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) append

asdoc reg mmxx l.lxx i.year if emajrelnewtmus==1, cluster(yccode_group) robust nested dec(3) tzok setstars(***@0.001, **@0.01, *@0.05) fs(10) font(CMU Serif) keep(l.wzsocdisx l.lxx _cons) stat(F) append
